Magnetocrystalline Anisotropy Energy of Transition Metal Thin Films: A 

Non-perturbative Theory 
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The magnetocrystalline anisotropy energy E an is of free-standing monolayers and thin films of Fe 
and Ni is determined using two different semi-empirical schemes. Within a tight-binding calculation 
for the 3d bands alone, we analyze in detail the relation between bandstructure and E ani3 , treating 
spin-orbit coupling (SOC) non-pertubatively. We find important contributions to E a nis due to 
the lifting of band degeneracies near the Fermi level by SOC. The important role of degeneracies 
is supported by the calculation of the electron temperature dependence of the magnetocrystalline 
anisotropy energy, which decreases with the temperature increasing on a scale of several hundred 
K. In general, E an is scales with the square of the SOC constant A so . Including 4s bands and s- 
d hybridization, the combined interpolation scheme yields anisotropy energies that quantitatively 
agree well with experiments for Fe and Ni monolayers on Cu(001). Finally, the anisotropy energy 
is calculated for systems of up to 14 layers. Even after including s-bands and for multilayers, the 
importance of degeneracies persists. Considering a fixed fct-Fe structure, we find a reorientation 
of the magnetization from perpendicular to in-plane at about 4 layers. For Ni, we find the correct 
in-plane easy-axis for the monolayer. However, since the anisotropy energy remains nearly constant, 
we do not find the experimentally observed reorientation. 
PACS: 75.30.Gw, 75.70.Ak, 75.70.-i, 73.20.Dx. 



I. INTRODUCTION 

The dependence of the total energy of a ferromagnetic 
crystal on the direction of magnetization originates from 
the magnetic dipole-dipole interaction as well as from 
spin-orbit coupling (SOC), as proposed by van VleckEl. 
The magnetic anisotropy energy is expected to be en- 
larged in systems of low symmetry, i.e. at surfaces, inter- 
faces, and thin filmsB or in one-dimensional systems such 
as quantum corralsu. Recently, a magnetization easy-axis 
perpendicular to the film plane has been observed for a 
wide variety of thin film .systems, for example for thin 
films of fee Fe on Cu(001)ca Some of these systems are 
promising candidates for magnetic high-density storage 
media. _ _ 

In spite of many theoretical attemptsQiL^I, the relation 
between the electronic structure and magnetocrystalline 
anisotropy energy E an i S could not be fully clarified so far. 
Some very important questions are subject to intense dis- 
cussion: (i) Which bandstructure details lead to signif- 
icant contributions to E an -i S ? Especially the treatment 
of degenerate bands, near the Fermi level has brought 
up controversiesEflJ. (ii) How does E anis depend on 
the SOC strength A so ? (iii) How is it influenced by the 
substrate lattice constant? Moreover, there is no unified 
thermodynamic and electronic theory to determine the 
temperature dependence of E an i S . Finally, the correct 
prediction of magnetic anisotropy for real systems still 
remains a challenge, since due to the quenching of or- 
bital angular momentum in 3d transition metal systems, 
Eanis is several orders of magnitude smaller than other 
contributions to the total energy of a crystal (typically 
about 0.1 — 1 meV per atom in ultrathin films). 

The magnetic anisotropy of thin films has been in- 



vestigated using two essenti^^different approaches. 
In semiempirical calculationsLji§E£Miil, the magnetocrys- 
talline anisotropy energy E an i S is determined by means 
of parametrized tight-binding bandstructures. Usually, 
spin-orbit coupling is restricted to second order pertur- 
bation theory. Qa_tlie|-Ojther hand, ab-initio calculations 
have been madeoOE3'E3 and lead to realistic bandstruc- 
tures. AlLcalculations make use of the controversial force 
theoremtj. Convergence, however, is difficult to achieve; 
sometimes, additional assumptions are made in-order to 
obtain converged results (state tracking methodti) . 

The structure of thin Fe films deposited on Cu(001) 
has been widely investigated, especially the dependence 
of the structure and magnetization orientation on the 
temperature. For films of less than 5 monolayers (ML) 
deposited at low temperatures, a distorted fec-structure 
is found, with magnetization perpendicular to the film 
plane. At 5 ML, a transition to in-plane magnetization 
is observed, as well as a restructuration of the film. It is 
still not clear if this reorientation transition is an effect 
of thp-stmctural changes taking place in the film at 4-5 

mimZtH 



In this paper, we investigate a simple quadratic Fe 
and Ni monolayer and fee multilayers systems up to 14 
ML epitaxially grown on the Cu(001) surface and neglect 
further interactions with the substrate. The bandstruc- 
tures are calculated within two different semi-empirical 
schemes, including SOC completely non-perturbatively 
without resorting to degenerate or non-degenerate per- 
turbation theory of any order. A tight-binding calcula- 
tion of the 3cZ-bands allows for a detailed, k-space re- 
solved analysis of the role of degeneracies for E an i S . It is 
shown that degeneracies located near the Fermi level can 
yield significant contributions, if they occur along lines 
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in k-space. We find for these that generally E an i S oc A 2 C 
holds. Including 4s-bands by means of the combined in- 
terpolation schemed and fitting the parameters to ab- 
initio calculations, we obtain the correct sign and val- 
ues of E an i s for the systems considered with this fully 
convergent method. That could be achieved neither by 
a fit using bulk parameters nor by employing a real- 
space density of states calculation, the so-called recursion 
methodO. Moreover, we find the characteristic scale for 
the temperature dependence of the magnetic anisotropy 
to be A so , rather than the bandwidth. This supports 
the significance of the lifting of degeneracies at Ep by 
A so and demonstrates the importance of contributions to 
magnetic anisotropy due to Fermi-edge smearing. 

Finally, we calculate the anisotropy energy of multi- 
layer systems. For systems of tetragonally-distorted Fe 
of 2 to 14 ML, we find a transition from magnetization 
perpendicular to the plane to in-plane magnetization at 
about 4 ML. We conclude from our calculation that the 
experimentally observed reorientation at 5 layers is not 
necessarily caused by a structural phase transition. For 
Ni, we find a nearly constant anisotropy energy from the 
fourth layer on, in disagreement with the results of Schulz 
and BaberschkeE3, who find a reorientation from in-plane 
to parallel magnetization at 7 ML. In both cases, the de- 
generacies near the Fermi-level are found to play an im- 
portant role for the dependence of the anisotropy energy 
on the film thickness. 

This paper is organized as follows: In Section II, the 
interpolation schemes (II. A, II. B, II. C) and the determi- 
nation of E an i S (II. D) are presented. The results for the 
tight-binding scheme for d-bands alone are shown in sec- 
tion III. A, the role of degeneracies is analyzed in detail 
in III.B while the results for the complete s- and d-band 
calculation for Fe and Ni monolayers on Cu(OOl) and 
other substrates are given in section III.C. The influence 
of crystal field splitting is investigated. Some aspects of 
the temperature dependence of magnetic anisotropy are 
considered in III.D, and the results for multilayer systems 
are presented in III.E and III.F. Section IV sums up the 
most important results. 



II. THEORY 

A. Bandstructures 

The magnetocrystalline anisotropy energy E an i S de- 
pends sensitively on the electronic structure of the sys- 
tem. To simplify the analysis, the bandstructure of the 
monolayer is calculated in two steps. First, the 3d-bands 
are described within a tight-binding scheme. Although 
the resulting E an i S as a function of the 3d-bandfilling n& 
shows already the most important features, the 4s-bands 
and s-d-hybridization have to be taken into account for 
a correct numerical evaluation of E an i S . 



For the 3d-bandS|, the tight-binding formalism intro- 
duced by FletchcrEil and Slater and KosterEa is adapted 
to the monolayer. The Hamiltonian H d = H at + AU is 
set up as a 10 x 10 matrix with respect to the basis of 
Bloch wave functions 

^ k (r) = -^^ e lk - R 0„(r-R). (1) 

R, 

Here, H at is the atomic Hamiltonian, AU the addi- 
tional crystal field in the monolayer, (pi, i — 1,...,5 
(i = 6, 10) are the atomic 3c? orbitals commonly de- 
noted by xy, yz, zx, x 2 — y 2 and 3z 2 — r 2 respectively, 
together with the spin eigenstate \\) (|J,)) with respect to 
the spin quantization axis zm- In the simple quadratic 
monolayer, only orbitals located on neighboring atoms 
are included. The extension to second nearest neighbors 
does not lead to further insightEj. With the x- and y- 
axes oriented along axes connecting nearest neighbors in 
the monolayer, the spin-polarized Hamilton matrix has 
(within the three-center approximation) the form 





= £o4 


- A^j + 2£i(cos2£ + cos 277) - J ex /2 




Tjd 


= E -\ 


- 2B 2 cos 2£ + 2B 3 cos 2r] - j' ex /2 




Tjd 

H 33 


= E -\ 


- 2B 3 cos 2£ + 2B 2 cos 2rj - J' ex /2 




Tjd 

-H44 


= E -\ 


- A^ e + 2B 4 (cos 2£ + cos 2rj) - J' ex /2 


(2) 


Tjd 


= E -\ 


- + A^ c + 25 5 (cos 2£ + cos 2?y) - J el 


:/2 


Tjd 


Tjd 


= H$ 10 = Hf 9 = 2J3 6 (cos2£ - cos2r/) 





and 



Hfi = Hf_ 5 i _ 5 + J ex for i = 6, 10, 
Hfi = Hf_ 5>i _ 5 + J' ex for i = 7, 8, 9. 

Here, £ = ^k x a and 77 = \k y a are the normalized com- 
ponents of the crystal momentum k, a is the lattice con- 
stant of the simple quadratic monolayer. For qualitative 
results it is sufficient to use bulk values for the param- 
eters of the paramagnetic bandstructure Bi, the crystal 
field parameter Ap g y Ni , and the spin splitting parameters 

J ex and J' ex . For NL t he parameters are taken from. WeL - 
ing and CallawayHaO, for Fe from Pustogowa et alt3'EB. 
The Bi and A v are listed in the first column of table Q. 
We have used J ex —0.1 eV and J' ex = 0.4 eV for Ni and 
Jex = J' ex — 1-78 eV for Fe. Due to the higher symmetry 
in fee or bec bulk crystals, only one crystal field parame- 
ter Ap e (AKj) appears in the corresponding Fe (Ni) bulk 
Hamiltonian. For the monolayer, one would have to con- 
sider three different A because of the reduced symmetry, 
but these parameters are not known. Hence, only Ap e ^ Ni 

has been considered in Eq. (J2j> . The influence of further 
crystal field effects-on E an i S in the monolayer, which was 
stressed by BrunailJ, is investigated in section III.C. 

For a quantitative comparison with experiment, how- 
ever, 4s-states have to be included_(within the so-called 
"combined interpolation scheme"tll) due to the strong 
overlap and hybridization between 3d- and 4s-bands in 
3d transition metals. According to the pseudopotential 
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method by HarrisonES, the 4s-electrons are described by 
a set of plane waves 



^K 3 k(r) 



1 



i(k-K 3 )T 



where the Kj are a set of reciprocal lattice vectors. They 
have to be chosen such that at least the lowest eigenstates 
in the considered part of the two dimensional Brillouin 
zone (irreducible part, see below) are described. For sim- 
ple quadratic monolayers, this yields Ki = (0,0), K2 = 
f (1,0), K 3 = f (0,1), K 4 = 2£ (1,1), K 5 = 2£ (-1,0) 
and Kg = — (1, — 1). To maintain the symmetry of the 
problem (and thus the correct occurrence of band degen- 
eracies that turn out— to be very important for E an i S ), 
symmetry factors i^Eil have to be introduced into the 
Hamilton matrix. This leads to 



for i 
else 



HI, = (V>K s k|i?|^K 3 k) 

V oa +a(k- Ki) 2 

\K i K, I' 1 1' 1 

Voo, Via, V11 V12 and V02 are the Fourier components of 
the pseudopotential, a is the dispersion of the 4s-band. 
The symmetry factors are: 



Fi 


= 1 




F 2 


= sin 2£ 




F, 


J sin2?7 


for r) > 0, 


= l 


else 


Fi 


= F2F3 




F 5 


J sin 2rj 

{ 


for r] < 0, 
else 


F 6 


= F2F5 





The s-d hybridization H sd between states ^af parallel 
spins is calculated according to Hodges et alxA with the 
parameters B\ and B 2 - To obtain accurate parameters, 
we perform a fit to the full-potential linear muffin-tin or- 
bitals (LMTO) calculation, for a free-standing Fe mono- 
layer by Pustogowa et al£3 and to the linear augmented 
plane wave (LAPW) calculation for a Ni monolayer by 
Jepsen et alia. The resulting parameters are listed in 
table |l]. In order to to reduce the number of free pa- 
rameters in the fit, the d-band parameters Bi and A v 
are still taken from the corresponding bulk crystals (see 
above). To obtain correct d-band widths, however, the 
Bi are scaled with the fitted parameters and for 
the spin-up and spin-down bands, respectively. Finally, 
the s- and d-bandwidths and s-d-hybridizalion parame- 
ters are scaled with t according to HarrisonO to take into 
account the Cu surface lattice constant a: 



(3) 



with clq the surface lattice constant of Fe or Ni, to the 
corresponding hopping parameters, and q being -5 for 



the dd parameters, -2 for the ss parameters and -7/2 for 
the sd parameters. The in-plane lattice constant is taken 
to be that of the Cu-substrate for all considered systems 
(a = 2.56 rA). This is correct for Ni, which is known to 
have a large pseudomorphic growth rangeca. For Fe how- 
ever, both an in-plane nearest-peigiibor distance similar 
to that of Cu and a smaller oneelrc-a have been reported. 



B. Spin-orbit coupling 

Spin-orbit coupling (SOC) between the d-states, lead- 
ing to magnetocrystalline anisotropy, is introduced in the 
usual form as H so = \ so \ ■ s. It can be expressed by 
the components of the orbital momentum operator 1 in 
the rotated frame (xm, Vm, ZmP- Here, zm is the spin 
quantization axis, "vdiich is parallel to the direction of 
magnetization (9, 4>)lB. 



SO SO 

so so 



lx M + il 



I — il 

L X M "i/M 



(4) 



Expressed in the basis of Eq. (|1|) , H so is a matrix function 
of the magnetization direction (6,<p). The SOC constant 
X so is taken from the corresponding atom: X so — 70 meV 
for Ni and 50 meV for F<£3. 

Unlike in usual tight-bindintt^alculationsErEllljO, SOC 
is included non-perturbativelyEZI in our treatment. Thus, 
we obtain important new information on how E an i S scales 
with the SOC constant A so , which contributes to our 
analysis of the origin of E an i S in terms of bandstructure 
properties (see below). 



C. Multilayers 



We build up the Hamiltonian of a system of I layers 
by coupling I monolayer Hamiltonians H^^ no l together. 
The coupling of the layers is described within the tight- 
binding nearest-neighbor formalism used for the mono- 
layer. Because of the missing periodicity in z-direction, 
we obtain terms that depend only on £ and rj. For the 
sake of simplicity, we take only cr-bonds into account and 
obtain the following terms for the coupling of the orbital 
j of the monolayer i with the orbital k of the monolayer 



-1, H 



i,i+l 
3,k 



H, 
H. 



M+i 
22 

i,i+l 
33 



= —2B\ cos 2r/ 
= -2i?icos2£ 



with i=l..J. The (18Z x 18Z) coupling matrix thus has 
only elements in the (I — 1) (18 x 18)-blocks just above 
and below the diagonal. The parameter B\ is the same as 
used for the monolayers, but it yet has to be scaled to the 
interlayer distance of the tetragonally-distorted system, 
according to Eq. ( ||). We consider equidistant layers. 
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For Ni, we take into account the reported compression 
of 3.2% to scale the interlayer hoppingsEl. For Fe, we 
assume, an expansion of about 5% as reported by Miiller 
et al!B. 



D. Anisotropy energy 



The magnetic anisotropy energy per atom is defined as 



E n 



,(n) :=E tot (9 = 0;n)-E to t(e = 7r/2,<fo;n), (5) 



where E to t(9, </>; n) is the ground-state energy per atom 
with a total of n 3d- and 4s-electrons per atom, and the 
magnetization direction is denoted by (8, tfijt^l. The in- 
plane angle </>o is chosen such that the resulting |-E a ms| 
is the largest possible. At first, the anisotropic dipolc- 
dipole interaction is neglected, since it does hardly de- 
pend on the electronic structure. Nevertheless, it may 
be of the same order of magnitude as the magnetocrys- 
talline anisotropy resulting from SOC and will thus be 
included later to obtain quantitative results. The total 
energy per atom E to t (with the k-space resolved energy 
Eis) is given by 

E tot {6,ct>;n) = 1 £ E k (9, <j>; n) 

k 

= 7j E e »*( > W° ( E mic(9, 4>) - E F (9, 4>\n)). (6) 



m , k 



with N the number of atoms. fo(AE) is the Fermi- 
function at zero temperature and Ep(9,(j)]n) is the 
Fermi-energy which, for a given bandfilling n, is deter- 
mined self-consistently by 



-E F (9, fan)). 



m,k 



E m k(9, <fi) is the m-th eigenvalue with crystal momentum 
k and magnetization along (9, <fi) of the Hamiltonian 



for the monolayer in the tight-binding scheme and 

77 77S | Tjd I 77Sd I / / 

tlmono — tl + tl + H + tl so 

for the monolayer in the combined interpolation scheme. 
For multilayer systems, we have the following Hamilto- 
nian: 



H = HL. 



H' 



H, 



coupling 



In Eq. ( we use the so-called force theorem, the 
validity of which has been assumed in all calculations of 
the magnetocrystalline anisotropy so far. 

The complete Brillouin zone (BZ) summation over k is 
performed as a weighted summation over the irreducible 



part of the BZ (for an arbitrary direction of magneti- 
zation). For the (i-electrons with SOC, that means a 
summation over 1/4 of the BZ. About 2000 points of 
the 1/4 BZ are then sufficient to achieve convergence. 
Note that we do not have to exclude any ptarts of the 
BZ to obtain convergence, unlike Wang et a/E3. Adding 
s-electrons and s-d hybridization implies a coupling of 
non-SOC coupled states with the SOC-coupled d-states 
and results in a reduced symmetry. It is then necessary to 
perform the summation over 1/2 of the BZ. We then need 
150 000 points to obtain the correct fourfold symmetry of 
the in-plane anisotropy energy as a function of the mag- 
netization direction in the plane (cos40). Fortunately, 
the out-of-plane anisotropy energy E an i S as defined by 
Eq.( H), which is larger by two orders of m agn itude in 
our calculation {E™~f ane ~ 1.2 fieV for Fe)E3 already 
converges for about 7000 points, so that calculations for 
systems of up to 14 layers are feasible. 



III. RESULTS AND DISCUSSION 

A. Monolayers within the tight-binding scheme 

In Figs. [I] and ||, results for E an i S as a function of 
the 3<i-bandfilling rid are presented (solid lines) for the 
parameters of Fe and Ni monolayers, respectively. We 
use the lattice constant of 2.56 rA to simulate epitaxial 
growth on Cu(001). These figures demonstrate the cor- 
respondence between electronic structure and magnetic 
anisotropy and show that our method will yield conver- 
gent results for the whole transition metal series and for 
large (Fe) and small (Ni) exchange coupling. They will 
be analyzed in the following. Yet, the numerical value 
of E an i s for Fe and Ni monolayers cannot be extracted 
from these figures until the 4s-electrons are included (see 
III.C), since the exact 3d-bandfilling of the monolayers is 
not known. 

Splitting the spin-orbit coupling matrix H so into two 
parts, one of them (Hf^ r ) containing only coupling be- 
tween states of parallel spin, the other one (H^ Upar ) 
between states of opposite spin, and recalculating E an i S 
as a function of rid with either of the two matrices in- 
stead of H so itself, we obtain the curves E^ r is (rid) and 
Eantl Par ( n d) j respectively (Figs, [l] and||, dashed and dot- 
ted lines respectively). Note that to a good approxi- 
mation ElZ s ( n d) + Knl7 ar ( n d) ~ E ams {n d ) is valid. 
For Fe parameters, E a ^^ aT {rid) is very small due to 
the large exchange splitting J ex that completely sepa- 
rates the spin subbands. Thus, E'^ l t i l f ar {ri d ) is ineffec- 
tive and may therefore be neglected for further analy- 
sis. The curve E^ is (rid) ~ E anis (n d ) consists of two 
parts of equal shape, viz. for n d e [0; 5] (spin-up band) 
and n d € [5; 10] (spin-down band). In the case of Ni, 
E anls( n d) an d E^ll par '(rid) are of the same order of mag- 
nitude, since there is a considerable overlap between the 
spin-up and spin-down subbands. 
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The curves E an i s (nd) show a number of pronounced 
peaks (A, B, C, E, F in Figs. and H), the origin of which 
has to be clarified. Two possible,contributions to E an i S 
are discussed in the literatureEiTllj: (i) The SOC-induced 
shifting of occupied, nondegenerate bands leads to con- 
tributions to E an i s in second-order perturbation theory 
with respect to the SOC constant A so : E an is oc AL.. The 
first order vanishes due to time reversal symmetryo. (ii) 
The contribution of the lifting of degenerate bands, which 
are shifted linearly with A so , depends on the fraction of 
states in k-space influenced by the degeneracy. Whether 
this fraction is of the order of A^ , which would yieldt2l 
Eanis oc X^ , or this fraction is of lower order AH|d-,thus 
would yield important contributions to E an j S J 11 !! 13 ! has 
been a controversial question. Anyway, the scaling of 
Eanis with X so can present important information about 
the dominant contributions to E an is- Thus, it is very 
useful not to restrict calculations to second order pertur- 
bation theory as has been frequently doneEJEJ. Remark- 
ably, we find E an i S (nd) oc Aj 0| ibr most of the rid values 
in agreement with Wang et alx3- Unlike stated by those 
authors, however, this does not rule out contributions to 
Eanis of the lifting of degeneracies (ii). In section III.B 
and Fig. |^, we show explicitly that such contributions 
play a very important role for E an is in the monolayers 
considered. This is true as well for the multilayers (s. 
Figs. [l3[ |TJ, and [TJ| and the discussion in sections III.E 
and III.F). 

The dependence of E an i S on scaling of all d-electron 
hopping parameters with a common parameter t was 
checked. We found that the overall shape of the curves 
E a nis{ n ) w ill not change if t is varied. |-E a ms| increases 
for decreasing t (decreasing bandwidth). This leads to 
the general trend: | E an is | increases with increasing lat- 
tice constant a of the mcmplayer, since t is proportional 
to a" 5 (see section III.G)E[ 

B. The electronic origin of E an is 

In this chapter we discuss in detail how the magne- 
tocrystalline anisotropy energy can be related to the elec- 
tronic bandstructure. A 3<i-band degeneracy can make 
large contributions to E an is, if (i) it is lifted by SOC for 
one direction of magnetization (rj^- ) and remains for an- 
other (z* 4 ); (ii) it is located near the Fermi level Ep; 
(iii) it runs along a line in k-space and (iv) the degen- 
erate bands have no or very little dispersion along this 
line. Before showing that such degeneracies indeed oc- 
cur in the bandstructures, we estimate their contribution 
within a linearized 'bandstructure' (see Fig. |J). If Ep 
is situated below or above the two subbands, no con- 
tribution to Eanis results: AEanis = 0. The maximal 
contribution occurs when the degeneracy lies exactly at 
the Fermi level Ep and amounts to: 



since the fraction F of involved states in the irreducible 

quarter of the BZ is F = = 2X so The 

preferred direction of magnetization is zf.j. 

Thus, AEanis is proportional to A^ for a degeneracy 
that occurs along a line with the involved bands being 
non-dispersive along that line. This agrees with the scal- 
ing of Eanis observed above. In their (estimate of the 
contribution of degeneracies, Wang et al£3 implicitly as- 
sume that the degenerate bands are dispersive in either 
dimension of k-space. This would lead to F oc \ 2 so and 
Eanis oc Af and justify the exclusion of degeneracies 
from their calculation in order to improve convergence. 
In the light of our results, however, this assumption is 
incorrect and it neglects very important contributions to 
E 

^ams • 

In Fig. ^, some degeneracies are shown in the band- 
structure of the Fe monolayer. For example, the degen- 
eracy A that occurs for M || z and is lifted for M || x 
is located at the Fermi level for rid — 7, 6 (dotted lines 
in Fig. H|) and leads to the peak A in Fig. [|. It runs 
along a line in k-space, which is shown Fig. g. According 
to Eq. (0), with = 0.6 eV/f (taken from the band- 
structure), this contribution should be AE anis w 4meV, 
which agrees in the order of magnitude with the calcu- 
lated value E an is(nd — 7.6) = 6 eV. 

Several tests have been made to support that hypoth- 
esis. Excluding the states influenced by the degeneracy 
A (4.3% of the total of 3<i-states) from the calculation 
of Eanis, the height of peak A is reduced to 40%. The 
k-space resolved analysis of E an is{nd = 7.6) (Fig. ||) also 
shows clearly that E an is results from the states near the 
degeneracy. 

Analogous degeneracies are found in the Ni bandstruc- 
ture contributing to the peaks E and F (Fig. |^). Note 
that the lifting of degeneracies can favor in-plane as well 
as perpendicular magnetization. This is in contradiction 
to the, results of Daalderop et at for a Co(lll) mono- 
layeiO, who state that degeneracies should always favor 
perpendicular magnetization. 

Since the 3d-band degeneracies are so important for 
Eanis, we analyze in the following the occurrence and 
lifting of degeneracies in the bandstructure. It can be 
shown that, in terms of the basis of Eq. (|]J), the Hamil- 
ton matrix H d (Eq. (||)) has the simplest block diag- 
onal form with only four off-diagonal elements (ODEs) 
Hf 5 = Hii and, equivalently, Hg 10 = H? 9 . To find 
out which additional ODEs are introduced by SOC for a 
given direction of the magnetization M, we analyze the 
form of H so in Eq. (^). States with parallel spins are 
coupled if they contain equal orbital momenta with re- 
spect to the spin quantization axis zm, whereas states 
with opposite spins must show a difference of one in the 
orbital momenta to yield nonvanishing ODEs. The real 
space components of the atomic states </>j, i = 1, 5, are 
composed of eigenstates of l z with the eigenvalues (-2,2), 
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(-1,1), (-1,1), (-2,2) and 0, respectively. In terms of eigen- 
states of l x one has the eigenvalues (-1,1), (-2,2), (-1,1), (- 
2,0,2) and (-2,0,2), respectively. This yields a coupling for 
M || z within the groups of states ipi with i — 1,4, 5, 7, 8 
and with i — 2,3,6,9,10, and, in the case of M || x, 
within the groups of states ipi with i = 2,4,5,6,8 and 
i = 1, 3, 7, 9, 10, respectively. In both cases, the Hamilto- 
nian can be split into two 5x5 blocks, and subbands be- 
longing to different blocks will intersect. Between states 
of the same block, the degeneracies will usually be re- 
moved. Especially the subbands ipi and ipi (and, corre- 
spondingly, ip§ and ipr) change their roles if the magneti- 
zation is changed from z to x and vice versa, because the 
orbitals xy and yz have different orbital momenta with 
respect to the x- and z-axes. These subbands will thus 
be involved in the lifting of degeneracies by altering mag- 
netization and possibly, as shown above, yield important 
contributions to E an i s . In the case of Fe parameters, the 
situation is even simpler since coupling between states 
of opposite spin (ipi and tpj with i < 5 < j) can be ne- 
glected. 

As an example, peak A in the curve E an i S (n) of Fe at 
n = 7.6 (Fig. |l|) results from the degeneracy A (Fig. |4|) 
of the subbands corresponding to the states ip-j and 
(?/>9, 0io). Thus, it occurs for M || z, and is lifted for 
M || x, since in the second case the subbands belong to 
the same block of the Hamiltonian, whereas in the first 
they do not. 

As a conclusion, it has been shown that 3d-band degen- 
eracies along lines of constant energy result in important 
contributions to E an is if they occur near the Fermi level. 
They can favor in-plane and perpendicular magnetization 
and need not occur near high symmetry points of the BZ. 
Thus, for (001) layers, it is not sufficient to consider only 
bands-at high symmetry points as was done by Daalderop 
et al£3 for a Co(lll) monolayer. Furthermore, for such 
contributions from degeneracies, E an i S oc X 2 so and, ap- 
proximately, E an i s oc 1/Jj£7 i s valid (the band dispersion 
is approximately proportional to the scaling t of the 
hopping parameters) which agrees with the observations 
reported above. Note that the analysis is very simple due 
to the analytic form and low dimension of the 3d tight- 
binding matrix, which is an advantage of the semiempir- 
ical scheme. It remains valid if the extension to s-states 
is performed (see below). 



C. The results of the combined interpolation scheme 

Results for E an i S (n) obtained from the combined in- 
terpolation scheme (including s- and d-bands as well as 
s-d hybridization) for the monolayer are presented in 
Fig. |g| for Fe parameters and Fig. [?] for Ni parameters 
with the lattice constant of the Cu(001) surface in both 
cases (solid curves; the discussion of the curves for two 
and three layers is postponed to section III.E and III.F). 
These results for the monolayer are similar to the curves 



for d-bands only (Figs [l] and ||) . n is the total filling of 
the s- and d-band (n = 8 for Fe and n = 10 for Ni). We 
find for a Fe monolayer i? a „i S (Fe/Cu) = —0.41 meV per 
atom, for Ni i? anis (Ni/Cu) = 0.10 meV per atom. The 
dipole-dipole interaction is included under the assump- 
tion of a point dipole located at each site, carrying the 
magnetic moment of the unit cell. The (spin) magnetic 
moment per atom is calculated from the bandstructure 
(m(Fe/Cu) = 3.3/i B and m(Ni/Cu) = 0.91^s)- The 
dipole anisotropy (equivalent to the shape anisotropy in 
the monolayer) always prefers in-plane magnetization. 
Alltogether, we obtain for the total magnetic anisotropy 
energy per atom of a Fe and Ni monolayer with the lattice 
constant of Cu(001) 

£*°„ t is (Fe/Cu) = -0.17meV and 
£*°™ s (Ni/Cu)= 0.12 meV 

with the easy axis perpendicular to the monolayer for 
Fe and in-plane for Ni. Note that corresponding ab- 
initio resulte-Jor a free-standing Fe-monolayer yielded 
—0.42 meV&U, but previous tight-binding calculations 
gave the too large value of —5.5 meVliil. 

In the case of Fe, the perpendicular easy axis of ultra- 
thin Fe films on Cu(001) is reproduced correctly. Direct 
comparison with a Fe monolayer on Cu(001) is difficult 
due to film growth problemso. It is common use to sepa- 
rate the anisotjpcy energy of thin films in a volume and 
a surface termaH 

2K 

E anis {d) = K v + — ^ (8) 
d 

The first term, K„, describes the thickness indepen- 
dent contributions to the anisotropy energy, and the 
second, K s , the thickness dependent contributions and 
the surface effects. Fowler i-and Barth measure the 
following anisotropy constantsO: K. v = 0.132 meV/atom 
and K s = 0.11 meV/atom for the distorted fcc-films at 
100 K. The value K V +2K S = 0.352 meV/atom is com- 
parable to our result. This result has been calculated 
with the measured anisotropy field using the bulk sat- 
uration magnetization of bcc-Fe. For.JNi, our result 
also agrees very well with experimentsBj which yields 
E ams (m/Cu) = 0.125 meV at 300 K. The anisotropy 
constants K s and K„ are temperature dependent. Mea- 
surements of the anisotropy constants as-,a function of 
the reduced temperature have been madeEa, but the cor- 
rect extrapolation to T=0 K is not known yet. While in 
experiment, the values of K v and K s have to be com- 
pared at the same reduced temperature because of the 
thickness dependence of T c , the theoretical values are for 
K and thus independent of the difference of absolute 
and reduced temperature. 

Note that in Fig. [?], the curve E an i S ln) for the Ni mono- 
layer (solid curve) has zeros near n = 10. Hence, the nu- 
merical result for Ni is not very stable and the excellent 
agreement with experiment should not be overempha- 
sized. Nevertheless, for Fe and Ni, the sign and the order 
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of magnitude of E an i S turn out to be remarkably stable 
upon parameter variations: Sign changes do not occur 
upon variation of the pseudopotential and s-d hybridiza- 
tion parameters by as much as 40%. Moreover, we find 
in agreement with Wang et a/.£3 a perpendicular easy 
axis also for Fe monolayers taking (001) surface lattice 
constants imposed by substrates such as Pd, Ag, and V 
(2.77 rA, 2.89 rA, and 3.03 rA), respectively. This sta- 
bility again demonstrates the validity of our results for 
Eanis- The good agreement of the results both with ab- 
initio theories and experiments is due to the fact that 
the parameters were obtained by a fit to ab-initio calcu- 
lations for Fe and Ni monolayers rather than taking bulk 
parameters. 

To investigate crystalr£eld effects, an additional pa- 
rameter A is introduced^ to take into account the dif- 
ferent effect of the monolayer geometry on orbitals that 
lie in the plane of the monolayer (xy and x 2 — y 2 ) and 
out-of-planc orbitals (yz, zx and 3z 2 — r 2 ). Additional to 
Eq. (H), the on-site energies of the latter are lowered by 
A with respect to the first. The dependence of E an i S on 
A is shown in Fig. || for Fe and Ni parameters (solid and 
dashed curve, respectively). Remarkably, A = 0.2 eV 
changes the sign of E an i S for both systems considered. 
Thus, it is important to determine A from the ab-initio 
bandstructures. In the case of Fe, the fit of the 3d-bands 
near the T-point of the BZ can be significantly improved 
by chosing A = 0.08 eV. The resulting E anis {F&/Cv) 
amounts to —0.30 meV, still with perpendicular easy axis 
even if the dipole-dipole interaction is added. For Ni, the 
introduction of A does not improve the fit. Those re- 
sults for A-.differ substantially from A = —0.5 eV given 
by BrunoEj which has been determined by a fit to the 
Ni(lll) monolayer but employed fat,both Fe and Ni(001) 
monolayers also. Pick and DreyssetU state that for (001)- 
monolayers a crystal field parameter is not necessary. For 
Ni, this is supported by our result; even in Fe, our value 
of A is small compared to other bandstructure parame- 
ters. Cinal et alH report A = -0.14 eV for the Ni(001) 
monolayer. 

Filially, a detailed investigation of the bandstruc- 
turesE3 shows that the analysis given in section III.B for 
3d-bands only is still valid for the combined interpola- 
tion scheme. As an evidence, consider Figs [j] and |6| (solid 
curves): There is a one-to-one correspondence between 
the peaks in E an i S in both curves. This correspondence 
can be shown to result from similar bandstructure details. 
In particular, the role of 3d-band degeneracies stressed in 
section III.B remains the same in the complete scheme. 



D. Temperature dependence 

One of the greatest challenges in the investigation of 
magnetic anisotropy is the calculation of reorientation 
transitions with temperature. Up to now, a complete 
electronic and thermodynamic theory is lacking. Here, 



one-particle effects of temperature are investigated. It 
turns out that they again support the role of degenera- 
cies for magnetic anisotropy and, moreover, are compara- 
ble in order of magnitude with the many-particle aspects 
usually considered^]. 

The free magnetic anisotropy energy F an i S depends on 
temperature T due to (i) the Fermi distribution of elec- 
tronic states /r(A£J), (ii) the hopping integrals, which 
depend on T because of the lattice expansion of the sub- 
strate, (iii) the entropy S(T) and (iv) the effects of spin- 
waves, resulting in a temperature dependence of the mag- 
netization M(T). In this work, the first three effects are 
analyzed. More precisely, the thermal expansion (ii) of 
the lattice constant a(T) is included by means of the em- 
pirical law a(T) = a(T = 0)(aT +1). a = 2 j-L0~ 5 /K is 
the expansion coefficient for the Cu substrata^. The ex- 
pression for the entropy (iii) of non-interacting particles 
is: 



S = 



k B ^ ( n "»k) hi (n mk ) + (1 - (n mk )) ln(l - (n mk )) 

771, k 



with (nmk) = f T (E mk (6, fa) - /i(6, fan)). In analogy 
to Eq. (0), the free magnetocrystalline anisotropy en- 
ergy F an i s is defined as the difference in the free energy 
F = E — TS for two different directions of magnetization. 

Fig. 10 shows F an i S (T) (<i-band calculation for the 
monolayer, Fe parameters, = 6). Including only 
Fermi statistics (i; dashed curve), the characteristic en- 
ergy scale for the decrease of |F „j s | with T is about 
1000 K (100 meV), which corresponds to the energy 
2 A so , but not to the 3d-bandwidth of approximately 3 eV. 
This becomes immediately plausible if one notices that 
the SOC-induced lifting of degeneracies occurs near the 
Fermi level. Thus, one expects a measurable effect on 
F an i s due to Fermi statistics as soon as fc^T becomes 
larger than or comparable to 2A so . In addition, we must 
conclude from our results that shifting of subbands far 
below the Fermi-level is not so important, since then 
Fanis could not be essentially lowered on such a small 
temperature scale. 

The characteristic increase of |-F a ms| with increasing 
temperature for T < 500 K is a direct result of the lifting 
of degeneracies. Consider again Fig. ||: for M || zf 4 
(lifted degeneracy), which is the energetically favored 
case, Fermi statistics induces only little changes in the 
occupation of the electronic states, if ksT < A so ; for the 



), however, states in the up- 



degenerate bands (M || z^ 
per band are significantly occupied even for ksT < X so . 
Thus, the total energy for M || z^ rises with respect 
to T = in this temperature range. This leads to an 
increase of |-F a nis| with increasing T, if fc^T < A so = 
50 meV (T < 500 K). 

The inclusion of lattice expansion (ii; solid curve in 
Fig. [lC]) has only a small effect on F an is ■ The narrowing 
of bands with increasing temperature due to the scaling 
of the hoppings leads to an increase in |-F tt nis| for small 
T, which was already discussed for T = 0. For larger 
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T, the influence of Fermi statistics on narrowed bands is 
larger, leading to a stronger decrease of |F arlis |. 

The entropy (iii; dotted curve) has a damping effect on 
the curve F an i a (T), but maintains the features discussed 
above. This results from the fact that, in the case of 
degenerate bands, the entropy is larger than for nondc- 
generate bands, since states located nearer to the Fermi 
level have larger entropy. 

FanisiT) was also calculated for bilayers (Fig. 11), 



taking into account all three mentioned effects and 
shows a decrease with increasing temperature on the 
same scale as for the monolayer. Hence, this anal- 
ysis of F an i S (T) shows the significant contribution of 
temperature-induced changes of the degeneracies to the 
anisotropy energy. It is remarkable that the three tem- 
perature effects mentioned above, and particularly the 
electron temperature dependence of the Fermi-function, 
are of equal magnitude as the temperature effects of spin- 
waves on M(T). 



E. Fe Multilayers 

Fig. [12] shows the calculated magnetic anisotropy en- 
ergy for Fe films of 1 to 14 layers. Calculations for both 
1/4 BZ and 1/2 BZ are included. The values obtained 
when the summation over k is performed over 1/4 BZ 
lead to periodically recurring positive values of E anis (for 
films of 2, 6, 9 and 12 layers). The positive value for film 
of 2 layers can be traced back to the occurence of degen- 
eracy A at the Fermi-level. For the other positive values, 
easy and hard axis are found to be in-plane, an effect of 
the wrong symmetry resulting from the summation over 
1/4 BZ, and the following overestimation of the in-plane 
anisotropy E™~f ione . A new degeneracy D is found to 
be responsible for the negative values. In Fig. we 
show the band structure of the Fe monolayer calculated 
within the combined interpolation scheme. The degen- 
eracy A observed for the monolayer in the tight-binding 
scheme is easy to recognize, and a new degeneracy D 
is found near the M-point for M || x. Degeneracy D is 
lifted for M || z, thus leading to a negative anisotropy en- 
ergy. The k-space resolved anisotropy energy shown in 
Fig. [l4| confirms the importance of degeneracy D, which 
causes the ring-shaped dip around M. The structure seen 
along the line LL' (s. Fig. ||) is the onset of the positive 
peak in the anisotropy energy caused by degeneracy A 
(s. Fig. |6|). Summation of the contributions of the k- 
points in the tenth of the BZ near M already gives half 
of the total anisotropy energy. Multilayer systems show 
per se more degeneracies than monolayers, and the con- 
tribution of these to the total anisotropy energy is not 
as clear as for the monolayer. Still, for a three layer 
system, we find again degeneracy D at the Fermi level, 
and recognize in the k-space resolved anisotropy energy 
(Fig. |l5|) the characteristic structure it causes around 
M. We perform the k-space summation again, this time 



over 1/2 BZ, thus respecting the symmetry of the s-d 
hybridized system. This reduces the importance of de- 
generacy A and we find positive values only at 2 and 6 
ML. Taking Fig. |l2] again and excluding the points of 
wrong symmetry (easy and hard axis in-plane) and the 
points where we find degeneracy A at the Fermi-level, 
we obtain the curve shown in Fig. ~E an i S is plotted 
as a function of I /I. We expected a linear behavior (s. 
Eq. ( ^)) and thus performed a linear least-square fit to 
the data. We then obtain K„ = -0.17 meV per atom and 
K s = -0.28 meV per atom, which is in very good agree- 
ment with Fowler und BartbH. 

Including the dipole-dipolc atusotropy energy as cal- 
culated by L. Szunyogh et altB, we would expect for 
Fe a change of the easy-axis from perpendicular to in- 
plane at 4 ML (#ff is (4 layers)=0.59meV). Our result 
indicates that the experimentally observed transition at 
5 ML might be an intrinsic quality of fct-films grown 
at low temperature. L. Szunyogh et al. calculated the 
anisotropy energy of thin fcc-Fe films on Au (001) and 
also observed oscillations. They obtained a reorientation 
transition from perpendicular to in-plane magnetization 
at 4 ML. 

The dependence of E an i S of Fe on the 3d- and 4s- 
bandfilling n is shown in Fig. 0. For the monolayer at 
a s- and d-bandfilling n=8, we are near a zero of the 
curve, and at n=8.2 we already have a positive value of 
Eanis caused by the growing influence of degeneracy A. 
We would thus expect a monolayer of a Fe x Co\_ x alloy 
to have an in-plane magnetization already at small Co 
concentrations. This was in fact measured by Dittschar 
et aEn for x=0.95. We would predict an increase of the 
anisotropy energy with increasing Co concentration. For 
3 layers, we would expect the same behavior, the struc- 
ture of the curve E a „i S (n) near n=8 being similar to that 
of the monolayer. This alloying behavior found both 
theoretically and experimentally supports the relevance 
of degeneracies for the anisotropy energy, as claimed by 
Daalderop et al. and disputed by Wang, Wu and Free- 
man. In this case, there is no doubt that the magnetic 
moment persists. 



F. Ni Multilayers 

The magnetic anisotropy of Ni calculated for systems 
of 1 to 14 ML is shown in Fig. [?]. We include again cal- 
culations using 1/4 of the BZ and 1/2 BZ, but this time 
no point has to be excluded. E a „i S of the second layer 
is much bigger than that of the monolayer, a fact which 
indicates that the influence of the substrate maybe can- 
not be neglected. The anisotropy then sinks again and 
remains approximately constant at a value of about 0.14 
meV (which is still bigger as the value, obtained for the 
monolayer). Schulz and BaberschkeE3 report for Ni a 
transition from in-plane to perpendicular magnetization 
at 7 ML, due to a large K„ which favors perpendicular 



8 



orientation of the magnetization. Our theory does not 
reproduce this reorientation. 

For Fe, the behavior of the films as a function of thick- 
ness could be related to the degeneracies occurring at 
the Fermi-level. The contribution of these degeneracies 
to the total anisotropy of the film would be expected to 
decrease with increasing number of layers, as their weight 
in the summation over all atoms (number of points in the 
BZ x number of layers) decreases: that is in fact what we 
find for Fe. For Ni, the contribution of the degeneracies 
to the anisotropy energy is not so evident, the minority 
and majority spin bands mix much more than in the case 
of Fe because of the small exchange coupling. This is a 
possible reason for the nearly constant anisotopy energy 
we obtain. The occurrence of a degeneracy of a /-times 
degenerated band would also probably lead to a thickness 
independent contribution to E a „i S . 

So far, no other monolayer calculation lead to the cor- 
rect in-plane anisotropy for the Niimonolayer. In a calcu- 
lation for the fct-bulk, ErikssoncS finds a perpendicular 
easy-axis, which is correct for fct-Ni, but wrong for fec- 
Ni. We obtain the correct in-plane anisotropy for the 
monolayer, but the wrong K„. So, a 3D calculation for 
fct-Ni does not really tackle the problem and explain the 
behavior of the magnetization. 

Ni is a delicate system. Maybe many-body effects can- 
not be neglected (i.e. the force theorem does not work 
well). However, a total-energy calculation made by Eriks- 
sont-H still yields the wrong sign for E an is for fcc-Ni bulk. 
The dependence between the anisotropy energy and the 
bandstructure seem to be very subtle and the smallest 
details can influence the results. 



IV. CONCLUSIONS 

A calculation of the magetocrystalline anisotropy en- 
ergy E an i S of Fe and Ni monolayers on Cu(OOl) is per- 
formed. In agreement with experiments, we find a per- 
pendicular easy axis for Fe and an in-plane easy axis for 
Ni. The results are fully converged without any addi- 
tional assumption to improve convergence. SOC is in- 
cluded non-perturbatively. It is an important result that 
large contributions to E anis can result from the SOC- 
induced lifting of degeneracies occuring along lines in 
k-space at the Fermi-level. The contributions of those 
degeneracies scale with the square of the SOC constant 
A so , as contributions from nondegenerate bands do. The 
occurrence and lifting of degeneracies in the 3d-band has 
been discussed in general. Evidence for the important 
contribution to E a „i s of the degeneracies at the Fermi- 
level are (i) the groove and the ring-shaped dip in the k- 
space resolved anisotropy for the monolayer in the tight- 
binding scheme and in the combined interpolation scheme 
respectively, (ii) the temperature dependence (the char- 
acteristic energy scale for the decrease of the free mag- 
netocrystalline anisotropy energy |F a „i s | as a function of 



the temperature is determined by A so ), (iii) the finite 
anisotropy energy at T c and (iv) the alloying behavior of 
Fe x Co\- x . We obtain for Fe a reorientation transition 
from perpendicular to in-plane magnetization at 4 ML, 
reorientation which is independent of any restructura- 
tion of the fct-film. Since it can be seen from Fig. [| and 
Fig. that both Fe and Ni do not exhaust the maximal 
anisotropy possible, our calculation of E an i S should also 
be important for the technologically relevant maximiza- 
tion of magnetic anisotropy by appropriate surface-alloy 
formation. 
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TABLE I. Bandstructure parameters within the combined 
interpolation scheme for Fe and Ni (OOl)-monolayers with 
lattice constant a. The parameters Bi and A v are taken 
from Pustogowa et al. 2S ' 29 for Fe and from Weling and 
Callaway 26 ' 27 for Ni (bulk parameters) . The other parameters 
are obtained from a fit to ab-initio calculations for freestand- 
ing (OOl)-monolayers by Pustogowa et al. 31 for Fe and Jepsen 
et al. 32 for Ni. 
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FIG. 1. Dependence of the magnetocrystalline anisotropy 
energy E a nis on the 3d-bandfilling rid for a monolayer with pa- 
rameters referring to Fe, calculated within the tight-binding 
scheme (solid curve). Negative values of Eanis yield perpen- 
dicular anisotropy. The origin of the peaks denoted by A, B 
and C can be traced back to degeneracies in the bandstructure 
(see text and Fig. The dashed and dotted curves show the 



contributions E^~ s and E^ll par to E anis from the spin-orbit 
coupling between parallel spins and antiparallel spins, respec- 
tively. 
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FIG. 2. Dependence of the magnetocrystalline anisotropy 
energy E anis on the 3d-bandfilling rid for a monolayer with pa- 
rameters referring to Ni, calculated within the tight-binding 
scheme (solid curve). Negative values of E an is yield perpen- 
dicular anisotropy. The origin of the peaks denoted by E and 
F can be traced back to degeneracies in the bandstructure 
(see text). The dashed and dotted curves show the contribu- 
tions E%£ u and E a a ^ ar to E anis from the spin-orbit coupling 
between parallel spins and antiparallel spins, respectively. 
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FIG. 3. Occurrence (full lines) and lifting (dashed lines) 
of a "line" degeneracy for two different directions of mag- 
netization and 2jjf, respectively, ki corresponds to one 
particular direction in k-space. Perpendicular to ki the in- 
tersecting bands are non- dispersive throughout the BZ. Note, 
the energy gained by the lifting of this degeneracy is given by 
AEanis = ^A so • F, if Ef falls in between the two subbands 
(dotted line). Here, F is the fraction of the involved states 
in k-space. Apparently, if Ef lies below or above the two 
subbands, AE a nis is zsro. 
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FIG. 4. Bandstructure of the 3d minority spin band of 
the Fe monolayer, calculated within the tight-binding scheme. 
The magnetization M is directed along the layer normal z 
(upper part) and in-plane along x (lower part). The degen- 
eracies denoted by A, B and C contribute to the peaks A, B 
and C in Fig. 1. The dotted lines denote the Fermi level for 
n d = 7, 6, respectively. T = (0, 0), X = (ir/a, 0), Y = (0, n/a) 
and M = (n/a, n/a) are the high symmetry points of the ir- 
reducible part (0 < k x ,k y < n/a) of the Brillouin zone, a is 
the lattice constant of the monolayer. 



FIG. 5. k-space resolved magnetocrystalline anisotropy en- 
ergy E anis (k,n d ) = E k (6 = O;n d )-£ k (0 = n/2,(j> = 0;n d ) in 
the irreducible part of the Brillouin zone (0 < k x ,k y < n/a; 
a is the lattice constant of the monolayer) for the bandstruc- 
ture of the Fe monolayer, n d = 7.6 electrons per atom in the 
d-band, calculated within the tight-binding scheme. Positive 
values of E a nis(k, n d ) favor in-plane magnetization. 
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FIG. 6. Magnetic anisotropy energy of Fe as a function the 
s- and d-bandfilling for 1 layer (solid curve), 2 layers (dashed) 
and 3 layers (dotted). Peaks A and D are caused by the 
respective degeneracies in the bandstructure shown in Fig. [l^. 
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FIG. 7. Magnetic anisotropy energy of Ni as a function the 
s- and d-bandfilling for 1 layer (solid curve), 2 layers (dashed) 
and 3 layers (dotted). 




a jL 

FIG. 8. Irreducible part of the two-dimensional Brillouin 
zone of Fe for the tight-binding scheme, a is the lattice con- 
stant of the monolayer. The main contribution to E an is at 
n — 8.8 (corresponding to = 7.6 in the tight-binding cal- 
culation) results from the lifting of degeneracies along the line 
LL'. 
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FIG. 9. Dependence of the magnetocrystalline anisotropy 
energy E anis on the crystal field splitting A for the Fe mono- 
layer on Cu(001), n = 8, (solid curve) and the Ni monolayer 
on Cu(001), n = 10, (dashed curve). Negative values of E an i a 
yield perpendicular anisotropy. The vertical line denotes the 
best fit for A for the Fe monlayer. In the case of Ni, the fit 
cannot be improved by the introduction of A (see text). 
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FIG. 10. Temperature dependence of F an is(T) for a 
Fe-parametrized d-band calculation for the monolayer with 
d-bandfilling n d = 6. For the dashed curve, only Fermi statis- 
tics is taken into account, for the solid curve the lattice ex- 
pansion is added, and the dotted curve includes the effects of 
Fermi statistics, lattice expansion and entropy. 
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FIG. 11. Temperature dependence of F an is(T) for 
a Fe-parametrized d-band calculation for 2 layers with 
d-bandfilling rid = 6. The calculation includes Fermi statis- 
tics, lattice expansion and entropy. 




Layers 

FIG. 12. Magnetic anisotropy energy of Fe as a func- 
tion of the number of layers calculated in the combined inter- 
polation scheme. The calculation for 1/4 BZ (dashed lines) 
yields periodic oscillations caused by the incorrect symmetry 
of E^ l 7 s p!a " e - Summation over 1/2 BZ corrects this problem. 
• and x are calculations with 15356 and 108228 points in the 
1/2 BZ respectively. 
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FIG. 13. Monolayer bandstructure of the 3s- and 4s-band 
for Fe-parameters, calculated within the combined interpola- 
tion scheme with the magnetization M parallel to the layer 
normal z in the upper part and in-plane parallel x in the lower 
part. High symmetry points are the same as in Fig. H. 
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FIG. 14. k-space resolved anisotropy energy for the Fe 
monolayer at a 3s- and 4d-bandfilling of 7.8. The ring-shaped 
dip near the M-point is caused by degeneracy D. 
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FIG. 15. k-space resolved anisotropy energy of Fe for 3 
layers at a 3s- and 4d-bandfilling of 8.0. The negative peak 
near the M-point is caused by degeneracy D. 



0.35 



0.3 



0.25 - 



0.2 



0.15 - 



0.1 




10 



12 



14 



Layers 



FIG. 17. Magnetic anisotropy energy of Ni as a function 
of the number of layers, calculated within the combined in- 
terpolation scheme. □ is a calculation for 1/4 BZ with 1722 
points, the solid curve is a calculation for 1/2 BZ with 6188 
points. 
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FIG. 16. Magnetic anisotropy energy of Fe as a func- 
tion of the l/l (I: number of layers). Including the 
dipole-dipole anisotropy energy, we obtain in-plane magneti- 
zation from the fourth layer on. Linear least square fit yields 
K„ = -0.17 meV per atom and K s = -0.28 meV per atom. 
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